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Abstract 

We use the Bohzmann equation in the relaxation time approximation to study the expansion of 
a dilute Fermi gas at unitarity. We focus, in particular, on the approach to the hydrodynamic limit. 
Our main finding are: i) In the regime that has been studied experimentally hydrodynamic effects 
beyond the Navier-Stokes approximation are small, ii) mean field corrections to the Boltzmann 
equation are not important, iii) experimental data imply that freezeout occurs very late, that means 
that the relaxation time remains smaller than the expansion time during the entire evolution of 
the system, iv) the experimental results also imply that the bulk viscosity is significantly smaller 
than the shear viscosity of the system. 
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I. INTRODUCTION 



The dilute Fermi gas at unitarity is a strongly correlated quantum fluid that serves as 
a new paradigm for strong correlations in many other systems, like dilute neutron matter, 
the quark-gluon plasma, and high Tc superconductors ^miM- The interest in the unitary 
Fermi gas derives in large part from the fact that the system provides a particularly simple 
realization of strong correlations. At unitarity the s-wave scattering length is inflnite, and 
the effective range and all other scattering parameters are zero. This implies that even 
though the system is very strongly correlated, details of the interaction are not important 
and the fluid is scale invariant. 

An important manifestation of strong correlations is the observation of nearly ideal hy- 
drodynamic flow Recently, a signiflcant amount of effort has been devoted to quantifying 
this observation by determining the shear viscosity t] or the dimensionless ratio of shear vis- 
cosity to entropy density r]/s |6|-[l0|. This effort is inspired, in part, by the conjecture that 
there is a universal lower bound on the shear viscosity to entropy density of strongly coupled 
fluids and by similar measurements of 77/ s for the quark gluon plasma created in heavy 



ion collisions at RHIC and the LHC 
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15|. 



Measurements of the shear viscosity of the unitary Fermi gas are based on using viscous 
hydrodynamics to analyze the expansion of a gas cloud after the conflning potential is 
removed. Alternatively, one can study collective oscillations of the cloud. The main difficulty 
with these analyses is that hydrodynamics breaks down in the dilute corona of the trapped 
gas, and that dissipative corrections from this regime can potentially be large. It is therefore 
important to study the crossover from weakly coUisional kinetic dynamics to almost ideal 
hydrodynamics. 

In this paper we use the Boltzmann (and Boltzmann-Vlasov) equation in the relaxation 
time approximation to study a number of issues related to the crossover from kinetics to 
hydrodynamics. In Sect. [Til we discuss the matching between the Boltzmann equation and 
the Navier-Stokes equation for a simple functional form of the relaxation time. We also 
study, for the system parameters realized in the recent experiments of Cao et al. 10|, the 
magnitude of hydrodynamic effects beyond the Navier-Stokes approximation. In Sect. IIIII 
we investigate the role of mean field effects in the Boltzmann-Vlasov equation. In Sect. IIVI 
we study the question whether there is any experimental evidence for freezeout, that is a 
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transition between hydrodynamic and ballistic expansion during the evolution of the system. 
In Sect. |V]we study a more realistic model of the relaxation time, and in Sect. I VI I we study 
the possible role of bulk viscosity. We should note that there is an extensive literature on 



using the Bo^ 
for example 



tzmann equation to describe the dynamics of trapped Bose and Fermi gases, see 
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24| . In this paper we focus, in particular, on the questions how close recent 



experiments are to the hydrodynamic limit, and how large the associated uncertainties in 
the shear viscosity are. 



II. BOLTZMANN EQUATION IN RELAXATION TIME APPROXIMATION 

We consider the limit in which the unitary Fermi gas can be described in terms of a 
quasi-particle distribution function f{x,v,t). The equation of motion for the distribution 
function is the Boltzmann-Vlasov equation 

^ + _ 1 djUmf + Ue.t) Of ^ ^ 

dt ^ dxi m dxi dvi 
where Uext = is the external confinement potential, Umf is a mean field potential, 

and C[f] is the collision term. We will study the effect of the mean field potential in Sect. IIIII 
For now we will set Umf = 0. Throughout this paper we will approximate the collision term 
using the relaxation time approximation 

C[f] ^ , (2) 

r 

where r is the relaxation time and fie is the local equilibrium distribution. The relaxation 
time does not have to be a constant, and we will study different scaling laws for r below. 
We are interested in the time evolution of a harmonically trapped cloud after the trapping 
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18| and use a scaling ansatz for the distribution 



potential is removed. We will follow 
function 

f{x,v,t)=T{t)fo{Rit),U{t)) (3) 

with 

R - ^ Ti - ^' ~ r - TT ^ (A\ 

where 6j, 9i are functions of t and /o(r, v) is a solution of the Boltzmann-Vlasov equation in 
equilibrium 

dfo _ 1 diUmf + Ue,t) dfo 



In the absence of a mean field potential this equation is solved by a Maxwell distribution. 
The scaling ansatz ([3]) breaks local thermal equilibrium only through the anisotropy of the 
temperature parameters 9i. The corresponding local equilibrium distribution fif, can be 
found by replacing 9i ^ 9 = (X^i 9i)/3. This distribution function is characterized by having 
the same mean kinetic energy as the non-equilibrium distribution /. In the free streaming 
limit r — !■ oo equ. ([3]) solves the Boltzmann equation exactly. In the presence of collisions we 
can obtain a differential equation for the parameters bi{t) and 9i{t) by computing moments 
of the Boltzmann equation. Integrating the Boltzmann equation over / d^U d^RUjRj (no 



sum over j) gives [17 1 



h + ^h-uf-^ = o. (6) 

The second term is due to the external potential and is not present if one considers an 
expanding cloud. Taking moments of the form / d^U d^R UjUj gives 

9, + 2^-19, = -I {9, - 9) . (7) 

Taking moments weighted with RjRj does not provide additional constraints. Together with 
the initial conditions bj{0) = 9j{0) = 1 and bj{0) = equations fl6|7|) describe the evolution 
of an expanding cloud. 

We expect the Boltzmann equation to reduce to hydrodynamics in the limit that the 
relaxation time is shorter than the expansion time, r ^ r^^p = {J2ibi/bi)~'^. In the limit 
r — 7- equ. ([7]) implies that 9i = 9 and 9 = (Hi The equation of motion for the scale 

parameter of an expanding cloud reduces to 

This equation is identical to the equation of motion that follows from the Euler equation in 
ideal hydrodynamics [s],!!^. If we keep terms of order r we recover dissipative terms in hydro- 
dynamics. For this purpose we compute the conserved charges (mass, momentum, and en- 
ergy) and the associated currents. The mass density is given by p = mn = m J d^vf{x, v, t). 
The momentum density is 

gi=m d^v Vif{x, V, t) = pVi (9) 



where Vi = {bi/b.i)Xi is the local fluid velocity. The energy density is 

1 /■ o .1 



S = -Jd^v mv^fix, V, t) = -pV^ + neo , (10) 



with eo = fT/e- The stress tensor is given by 

Uij = J mViVjf{x, V, t) = pViVj + 5ijp9i{v1) , (11) 

where {vf) = / d^v fo{R, v) with V = Jli bj is an average with respect to the equihbrium 
distribution in the local rest frame. We can write f = fie + Sf and split equ. f|TT]) into a 
local equilibrium and dissipative part. The dissipative part is 

SU,, = -d,,nTie - I) • (12) 

We can use equ. (JTj) to express 6i/6 in terms of 6i and 6j. The result is explicitly pro- 
portional to r, so that at leading order we can use the relation between 6i and 6j in ideal 
hydrodynamics. We find 

5n, = -5.,nT.r(2|-^E|) • (13) 

Using Vi = Xiipi/hi) we see that equ. ( fT3|) matches the viscous stress tensor in hydrodynamics, 
= -r]{diV, + djVi - pijid ■ V)) - (Sijid ■ V), if we set 

r] = nrTie , (14) 



and ( = 0. Equ. (IT^ is a standard result in kinetic theory 



22]. The bulk viscosity is 



expected to vanish in the dilute Fermi gas at unitarity |25|, |26|, but we will see that the 
relaxation time approximation leads to a vanishing bulk viscosity even in cases where we 
expect the physical result to be non-zero. Scale invariance also implies that rj = nan{y), 
where a„ is a function oi y = mT/n^^^, and r = Tn{y)/T. The matching condition (fl^ 
then implies that r„,(?/) = an{y)- The simplest case is that a„ is a constant. In this case the 
relaxation time is only a function of temperature r = a^/T/e = a„/(To^), where Tq is the 
initial temperature. 

In Fig. [1] we compare the solution of the Boltzmann equation to a solution of the Navier- 
Stokes equation. The Navier-Stokes solution is described in [2^, and the parameter a„ 
was fitted to data taken by Cao et al. 1^ at an initial energy E/Ep = 3.61, where Ep = 
(SNY^^uj. The experiment involves an axially symmetric trap with uj± = cOx = cOy- The trap 
frequencies are Uj_/{2'7i) = 5078.7 Hz and Uz/{2n) = 175.7 Hz, the number of particles is 

= 2.61 ■ 10^ and we have defined u = \^^^uj± with A = ujz/co±. The relaxation time was 
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FIG. 1: This figure shows the matching between kinetic theory and Navier-Stokes hydrodynamics. 
We show the evolution of the aspect ratio as a function of time. The sohd points are data 
taken at an initial energy E/Ep = 3.61 



10(1 . The solid line shows a solution of the Euler equation, 
the long dashed line is a solution of the Navier-Stokes equation where the viscosity coefficient 
On = 22.1 (r/ = a„?i) was adjusted to reproduce the data, and the short-dashed line is a solution 
to the Boltzmann equation in the relaxation time approximation with r = On/T. 

fixed according to the relation r„ = a„. Fig. [1] shows the time evolution of the aspect ratio 
^r(^) = b±(t)/bz(t). We observe that the difference between the kinetic theory calculation 
and Navier-Stokes hydrodynamics is indeed small, indicating that in the regime that has 
been studied experimentally higher order hydrodynamic effects are small. 

This issue is studied in more detail in Fig. O For a given value of the viscosity coefficient 
a„ we determine the Navier-Stokes evolution and then fit the relaxation time coefficient r„ 
in the kinetic calculation to provide the best fit to hydrodynamics. The figure shows the 
dimensionless parameter 

= {SNXy/^Eo/Ep ^^^^ 
that governs the hydrodynamic evolution as a function of the dimensionless relaxation time 

f. = ^. (16) 

In the hydrodynamic limit equ. f ll4p implies = t„/3. This equation is satisfied at small 
r„, but for t„ > 1 the effective shear viscosity is smaller than the one predicted by the linear 
approximation. The fit in Fig. [1] corresponds to f„ = 0.78, which is close to the regime 



0.5 1.0 1.5 2.0 

FIG. 2: This figure shows the relation between the relaxation time in kinetic theory and the shear 
viscosity. We plot the dimensionless relaxation time r against the dimensionless shear viscosity 
defined in the text. The points were obtained by first solving the Navier-Stokes equation with a 
given values of /3,, for a trapped cloud with initial aspect ratio Ar{{)) = 0.035, and then determining 
the value of r that leads to the best fit of the solution of the Boltzmann equation to the Navier- 
Stokes result for A{i{t) in the regime t < 1400 ms. The line show the leading order relation 
~ ^l"^- Note that the long dashed line in Fig. [T] corresponds to = 0.25. 

where higher order effects become large. Qualitatively, the behavior of as a function 
of f„ is easy to understand. At moderate r„ the relaxation time approximation implies 
that the dissipative stress tensor equ. f|T2|) does not reach the hydrodynamic limit equ. ( JT3|) 
instantaneously, but over a characteristic time r„/Tie. For an expanding system this means 
that the time integral of the dissipative stresses in the kinetic theory is smaller as compared 
to the hydrodynamic limit. At large r„ the solution of the Boltzmann equation approaches 
the ballistic limit, and dissipative effects saturate as f„ — ^ C)0. 

III. MEAN FIELD EFFECTS 

In this section we investigate possible mean field effects in the Boltzmann equation, in 
particular the question to what extent these effects cause uncertainties in the extraction of 
the relaxation time. Without mean field effects the equation of state is P = nT and the 
relation between energy density and pressure is £^ = |P. In the unitary gas scale invariance 
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FIG. 3: This figure shows the effect of the mean field potential on the evolution of the aspect ratio 
Aji{t). The solid line shows the result in ideal hydrodynamics, and the points are experimental 
data from |lo| taken at an initial energy E/Ep = 3.61. The short dashed line is the result of a 
relaxation time fit to the data without mean field effects. The band around the dashed line shows 
the variation of Aji{t) when mean field effects included. The upper and lower limits of the band 
correspond to ^ = ±0.1, respectively. 

restricts the equation of state to be of the form P = nTh(y) with y = mT/n^/'^ and h{y) 
is a universal function which approaches unity in the high temperature, low density limit. 
The relation £ = |P is exact irrespective of the functional form of h{y). Scale invariance 
restricts the form of the mean field potential to be Umf = g{y)n{x)°' with a = 2/3. For 
simplicity we will assume that g{y) ~ const. The mean field potential does not modify the 
equation for 9j. The equation for bj becomes 

s,..j.,-.;a + .|^g-^).o. (17) 

with V = Hj bj. The parameter ^ is defined by 

In the limit of ideal hydrodynamics (r — 0) we have 9j = 9 = V^^^^. This implies that 



^= i^.Jl7:nr (18) 



mean field effects do not modify the evolution a scale invariant, non-dissipative gas. This 

nn 

result is known from studies of the unitary gas in the context of ideal hydrodynamics |8l. I27l|. 
The basic observation is that the Euler equation only depends on the relation P{S), and 
not on the equation of state P{n,T). 
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For a finite relaxation time the relation 6j = V~^^^ is only approximately satisfied and 
mean field corrections play a role. In order to estimate the size of mean field effects we relate 
the parameter ^ to thermodynamic properties of the unitary gas. The static Boltzmann- 
Vlasov equation (|5]) imphes 

l^v') + ^Y^(f^w) = jT.-^nr^) ■ (19) 

For a scale invariant gas with a = 2/3 the LHS of this relation is equal to the internal energy 
of the gas. Equ. ( fT9l) is then the well-known Virial theorem 28l-l30|: The total internal energy 
of a harmonically trapped system is equal to the potential energy. We can use equ. ( |T9i) to 
express ^ in terms of the total energy of the gas. For a = 2/3 and in the limit ^ ^ 1 we find 

e = - 1 . (20) 

^ 3NT ^ ' 

For an ideal gas the total (internal plus potential) energy is Eq = 3NT and ^ = 0. We have 
computed ^ using the equation of state described in the appendix of |27j. This equation 
of state is a parameterization of the recent experimental results of Nascimbene et al. 31 |. 
We find that in the normal fluid regime |^| < 0.1, where ^ < at high temperature regime 
and ^ > near the critical temperature T^. We note that at unitarity corrections to the 
equation of state due to quantum statistics are parametrically of the same order as mean 
field corrections. At high temperature quantum corrections are numerically small, but near 
Tc equ. fl20|) is not reliable. 

In Fig. [3] we show the effect of mean field corrections in the regime —0.1 < ^ < 0.1 
on the evolution of the system. Positive values of ^ tend to push the evolution closer to 
the solution in ideal hydrodynamics, and therefore act like an effective negative viscosity. 
Negative values of ^ act an effective positive shear viscosity. Fig. [3] shows that mean field 
effects are small compared to dissipative effects. We can quantify this observation using the 
data at E/Ep = 3.61. A fit of the relaxation time using the Boltzmann equation without 



mean field effects gives f:„ = 1.14. For this energy the equation of state described in [27 1 
gives ^ = —0.014. If we refit the relaxation time using the Boltzmann- Vlasov equation with 
this value of ^ we find fn = 1.12, which is a ~ 2% correction. 
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FIG. 4: This figure shows the effect of freezeout on the time evolution of the aspect ratio. The data 
points show the results of Cao et al. also shown in Figs. [1] and [3l The solid line is the relaxation 
time fit shown in Fig. [3l The dashed lines show the same fit but with the scaling of the relaxation 
time changed from r ~ 1/T to r ~ T^/^V at a freezeout scale factor bfr = 10,5,1 (from top to 
bottom). 

IV. FREEZEOUT 

Up to this point we have considered a relaxation time of the form r = r„/T with r„ = 
const. In this case the relaxation time scales as r ~ 9^^ ~ V"^^^. The expansion time, on 
the other hand, scales as r ~ V^^^ (at very late time the expansion changes from being 
two-dimensional to three-dimensional and r ~ V^^'^). This implies that the system will 
eventually freeze out, and the nature of the expansion will change from hydrodynamic to 
ballistic. We note, however, that for the values of r„ that fit the data shown in Fig. [T|3l 
freezeout occurs late, for Af> > 1. At this time the density has dropped by a factor of order 

To maintain hydrodynamic behavior despite the large drop in density requires the system 
to be very close to the unitary limit. At unitarity the mean free path scales as Imfp ~ 
(ncr)~^ ~ VX']^^ ~ where XdeB ~ (mT)~^/^ is the thermal wave length. This implies 
that the the mean free path grows more slowly as the system size. For a weakly collisional 
system, on the other hand, we have Imfp ~ (na^)~^ ~ V and the mean free path grows more 
quickly than the system size. 
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In this section we wish to check whether the data support the idea that freezeout occurs 
late. For this purpose we impose a freezeout time tfr so that for t < tfr the relaxation 
time scales as r ~ t„/T and for t > tfr the relaxation time scales as r ~ 
is the behavior expected for a constant cross section a ~ a^. The parameter r„ is taken 
from the previous fit and 7 is adjusted so that r is continuous at t = tfr. In Fig. H] we 
show the evolution of the aspect ratio ^^(t) for different choices of the freezeout time. The 
solid line shows the result for tfr — )■ 00, and the dashed curves correspond to freezeout scale 
parameters b±{tfr) = 10, 5, 1. We observe that tfr — > 00 provides the best fit to the data, and 
b±{tfr) < 5 is clearly excluded. We conclude that for the conditions explored experimentally 
(A^ ~ 2-10^ and E/Ep < 3.61) the density can drop by a least a factor ~ 25 before freezeout 
occurs. 



V. HIGH TEMPERATURE LIMIT OF THE RELAXATION TIME 



{f{x,v,t)- fUx,v,t)) . (21) 
const the relaxation time 



In the high temperature limit the viscous relaxation time can be c omp uted reliably. The 
result is r = r]/{nT) with 77 = r]Q{mTf/'^ and r^o = 15/(320r) [221, M |32j . This implies that 
even within the relaxation time approximation the collision term is a nonlinear functional 
of the distribution function. We have 

C[f] = \jd'vf{x,v, 

We note that even in the simplest case 1/r ~ Tie/Tn with r, 
is a functional (via the temperature) of the distribution function. However, for the scaling 
ansatz in equ. ([3]) the temperature is only a function of time. The collision term in equ. (12T|) 
is more complicated because the relaxation time is a function of time and position (but 
not of velocity). This means that it is difficult to find exact solutions of the Boltzmann 
equation. We obtain an approximate scaling solution by using the same ansatz for f{x,v, t) 
as before. The parameters bi and 6i are again fixed by taking moments of the Boltzmann 
equation weighted with R^, Uf and RiUi. We find that only the U"} equation is modified. As 
a consequence equ. ([7]) is replaced by 



9, + 



(22) 



where (no) = / d^rnoir)"^/ J d^rnQ{r) is the average density in the initial state. Note that 
for a Gaussian distribution (no) = 2~'^/^no(0). We also observe that to leading order in r 
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the parameter 9 evolves as ^ = V~^/^. This imphes that the RHS of equ. f l22p scales as 
T(fi/Tn with r„ = r]o(mTo)^^'^ / (no) . As a consequence, the evolution of the system in the 
Navier-Stokes regime (that is to first order in r) is indistinguishable from the result shown 
in Fig. [U What is new is the explicit relation between r„ and the shear viscosity in the high 
temperature, low density limit. 



This relation can be compared to the analysis presented by Cao et al. in [10|]. In this 
work the authors attempted to extract r^o from experimental data on elliptic flow in the high 
temperature regime using the Navier-Stokes equation. The problem with the Navier-Stokes 
analysis is that a density independent shear viscosity leads to an infinitely large dissipative 
contribution from the dilute corona. Cao et al. argued that a correct treatment of the 
corona has to take into account the fact that at low density the dissipative stress tensor 
cannot instantaneously relax to the Navier-Stokes form. They proposed that this effect 
can be taken into account by using a simple model for the viscosity which is of the form 

more detailed justification for this model can be found in 
|8|. The model of Cao et al. can be written as t] = with = r]Q(mT)^^'^ / n{0) . Because 
the temperature drops (approximately) as T ~ n{0)^^^ this corresponds to a constant value 
of On = rio{mToy^^ /nQ(0). This equation for q;„ has the same structure as the result that 
follows from equ. ( l22l) . but it differs by the numerical factor (no) /no{0). This difference is 
clearly very significant, and it represents the most important uncertainty in the extraction 
of ?7o from data. This uncertainty can only be addressed by studying more accurate solution 
of the Boltzmann equation in the case of a density dependent relaxation time. 

VI. BULK VISCOSITY 

Finally, we wish to study the possible effects of bulk viscosity on the evolution of the 



cloud. At unitarity we expect the bulk viscosity to vanish [25|. There are nevertheless at 
least two reasons for studying dissipative effects associated with bulk motion. The first is 
to verify that the theoretical prediction ( = at unitarity is indeed correct. The second is 
to understand how bulk viscosity manifests itself as one moves away from unitarity. 

In Sect. |TT]we showed that the relaxation time approximation to the Boltzmann equation 
(without mean field effects) leads to a traceless dissipative stress tensor, see equ. f lT2|T3|) . We 
can break scale invariance by using a mean field potential with a ^ 2/3. However, within 
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FIG. 5: This figure shows the effect of shear and bulk viscosity on the evolution of the aspect ratio 
Aftit) in viscous hydrodynamics. The solid line shows the result in ideal hydrodynamics, the long 
dashed line corresponds to the Navier-Stokes equation with = 0.25, and the short dashed line 
shows the Navier-Stokes result for = 0.25. The system parameters are the same as in Fig. [TJ 

the approximations used in Sect. IIIIl scale breaking in the Boltzmann-Vlasov equation does 
not lead to a trace-term in the dissipative stress tensor. We will therefore study the role of 
bulk viscosity using the Navier-Stokes equation. We will follow 0, Q and use a scaling 
ansatz for the hydrodynamic variables. 



The functional form of the density n and the velocity field Vi agree with the results obtained 
in Sect. HIl In the Navier-Stokes limit the function aj(t) can be related to the parameter 6i{t) 
in the distribution function by Oj = Oiuf/hf. For an axially symmetric system the equations 
of motion for hi and are given by 






(23) 
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FIG. 6: This figure shows log(x^) contours for a hydrodynamic fit to the data of Cao et al. shown 
in Fig. [TJ The minimum of corresponds to the dark region. The two parameter and 
control the bulk and shear viscosity of the system. 

Here, (3ri is defined in equ. (fTSll and is the analogous parameter that controls the effect of 
bulk viscosity 

= {3NXy/^Eo/Ep ^^^^ 
with ( = C„n. For = and to first order in eqns. fl24ti27p are equivalent to eqns. fl6|71) . 
Solutions of these equations for f3(^) = (0.25, 0) and = (0, 0.25) are shown in 

Fig. |5l We see that the effects of shear and bulk viscosity are clearly distinguishable. Shear 
viscosity leads to a characteristic curvature in the time dependence of the aspect ratio. 
This is related to the fact that shear viscosity is more efficient in moving kinetic energy 
from transverse motion to longitudinal motion, combined with the fact that the longitudinal 
expansion takes place on a longer time scale. A comparison with the data of Cao et al. is 
shown in Fig. O The figure shows contours for a fit to the data at E/Ef = 3.61 using 
both and as fit parameters. We observe that the absolute minimum is at = 0.25 
and = 0.0 («„ = 22.1 and (n = 0.0). We conclude that the data at unitarity favor a 
vanishing bulk viscosity and a non-zero shear viscosity. 
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VII. SUMMARY AND CONCLUSIONS 



In this paper we have studied a number of issues related to the crossover between kinetic 
theory and hydrodynamics. In this section we briefly summarize the main findings and point 
to open problems. 

1. Dimensional analysis implies that the shear viscosity and the relaxation time can be 
written as rj = and r = Tn/T, respectively. Scale invariance dictates that at 
unitarity a„ and Tn are only functions oiy = mT jr?!"^ . The matching between kinetic 
theory in the relaxation time approximation and hydrodynamics is simplest in the case 
that and r„ are constant. In this case the matching condition is simply a„ = r„, 
and we showed that for systems that have been studied experimentally terms of higher 
order in r are not large. 

2. The Boltzmann equation describes the evolution of a system that obeys an ideal gas 
equation of state, P = nT. Corrections due to a more complicated equation of state 
can be taken into account using the Boltzmann- Vlasov equation and a mean field 
potential. At unitarity and in the limit r — )■ the mean field potential has no effect 
on the evolution of the system. As a consequence, mean field corrections remain small 
even if the relaxation time is not zero. Using the mean energy per particle from a 
realistic equation of state we find that the uncertainty in the shear viscosity related 
to mean field effects is ~ 2%. This is consistent with the effect due to the deviation 
of the equation of state from P = nT found in calculations using the Navier-Stokes 



equation [27 1. 



3. We also studied the effect of a possible freezeout, that means a transition from hy- 
drodynamic to ballistic behavior, on the evolution of the system. Freezeout could 
occur because the scattering length is not strictly infinite, or because non-equilibrium 
effects not captured by the Boltzmann equation cause the growth of the transport 
cross section to lag behind the equilibrium prediction. We showed that the experi- 
mental data obtained by Cao et al. [lO| require that freezeout has to occur late, for 
values of the scale parameters &± > 5. This implies that the unitary Fermi gas re- 
mains hydrodynamic despite the fact that the density drops significantly during the 
expansion. 
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4. Finally, we studied the possible role of bulk viscosity on the evolution of the system. 
We showed that there are characteristic differences between dissipative effects due to 
bulk and shear viscosity, and that the data are best described by shear viscosity only. 



There are a number of issues that remain to be resolved. We have studied the Boltzmann 
equation with a relaxation time that scales as r ~ T^^'^/n. This is the scaling which is 
predicted by the linearized form of the full collision terms using the scattering amplitude in 
the unitary limit. It corresponds to a shear viscosity that is only a function of temperature 
and not of density. For simplicity we have employed the scaling ansatz equ. ([3]) for the 
distribution function and solved for the coefficients by taking moments. We obtain an equa- 
tion of motion which is equivalent to the /T with an effective r„ = r]{TQ)/{no). 
This result is in qualitative, but not in quantitative, agreement with an estimate based on 
hydrodynamics with a finite relaxation time js, 10|. In order to resolve this disagreement we 
need to find more accurate solutions of the Boltzmann equation. For infinite systems there 
are solutions to the linearized Boltzmann equation that can be matched to second order 
hydrodynamics [32], but the corresponding solutions for an expanding system have not been 
determined. 

We also observed that within the approximations used in this work the Boltzmann equa- 
tion does not account for bulk viscosity even if scale invariance is broken by mean field 
corrections to the equation of state. It is known that the integral of the frequency depen- 
dent bulk viscosity is non-zero if the equation of state violates scale invariance 33|]. It is not 
know what the leading 1/a term in the bulk viscosity is, and what terms in the Boltzmann 
equation have to be kept in order to reproduce this term. 
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